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Statement  of  Problem 

^  The  problem  addressed  in  the  work  supported  was  to  develop  a  method  for 
accurate  all  electron  self  consistent  calculations  of  the  electronic  structure 
of  small  clusters  of  transition  metal  atoms  on  the  basis  of  local  spin  density 
functional  theory  and  using  a  basis  of  symmetrized  combinations  of  Gaussian 
orbitals.  Our  intention  was  to  use  the  clusters  to  study  the  effects  of  small 
concentrations  of  foreign  atoms  in  transition  metal  hosts  (example:  carbon  in 
iron),  because  calculations  for  clusters  could  be  made  self  consistent  to  an 
extent  not  normally  possible  for  impurities  in  bulk  solids.  Subsequently, 
experimental  observations  of  the  properties  of  iron  clusters  by  research 
groups  at  the  Exxon  laboratories  and  elsewhere  have  stimulated  us  to  consider 
free  clusters  in  more  detail. 


4 


Results 

We  begin  our  work  by  attempting  to  modify  the  standard  quantum  chemistry 
program  GAUSS76  to  perform  density  functional  calculations.  This  attempt  was 
successful  for  small  systems  (diatomic  molecules).  The  results  obtained  for 
LiH,  CO  and  NO  indicate  that  spin  density  functional  theory  leads  to  an 
excellent  account  of  the  ground  state  properties  of  simple  molecules. 

However,  it  did  not  prove  feasible  to  study  the  large  (lb  -  20  atom)  systems 
of  primary  interest  to  us  with  this  approach. 

Instead  we  developed  our  own  program.  The  essential  features  are  that  it 
makes  extensive  use  of  symmetry  to  enable  efficient  calculations  for  larye 
systems.  An  auxiliary  basis  was  introduced  for  variational  fitting  of  the 
charge  density  to  eliminate  the  necessity  for  the  computation  of  four  center- 
two  electron  integrals.  A  numerical  grid  was  developed  for  the  evaluation  of 
matrix  elements  of  the  exchange  correlation  potential.  We  believe  this 
program  is  the  most  accurate  available  at  the  present  time  for  cluster 
calculations. 

Results  were  obtained  for  the  following  clusters:  Feg,  Fejb,  Fe^C  (with 
body  centered  cubic  geometry)  and  Fe^,  Ni^,  and  Ni ^ 9  (with  face  centered 
cubic  geometry).  The  atomic  spacing  was  taken  to  that  of  the  bulk  metal  (the 
geometries  of  experimentally  observed  clusters  are  unknown).  Specific  results 
are  described  in  detail  in  the  publications  which  are  listed  subsequently.  In 
a  general  way,  they  may  be  summarized  as  follows: 

The  general  features  of  the  density  of  state  for  the  cluster  resemble 
those  of  the  bulk  metal  as  obtained  from  band  calculations  for  clusters  as 
small  as  15  atoms.  This  was  particularly  striking  for  the  case  of  the  Fe^ 
cluster:  Not  only  does  the  overall  width  in  energy  of  the  d  level  complex  in 


the  cluster  agree  with  that  of  the  bulk  but  the  positions  and  heights  of 
several  of  the  major  peaks  in  the  density  of  states  correspond  as  well.  Our 
results  for  clusters  are  in  much  better  agreement  with  those  for  the  solid 
than  are  those  of  other  workers  who  have  used  different  methods  (DVM,  MS-Xa), 
and  support  the  view  that  the  density  of  states  depends  only  weakly  on  the 
boundary  conditions  imposed. 

All  of  the  clusters  studied  were  ferromagnetic;  i.e.,  have  net  maynetic 
moments.  The  moments  are  larger  than  for  the  bulk,  roughly  3pg/atom  for  the 
iron  clusters.  This  has  turned  out  to  be  in  agreement  with  observation. 
However  the  spatial  distribution  of  the  spin  density  in  the  iron  clusters  is 
quite  different  from  that  of  the  bulk.  The  central  atom  has  more  minority 
(down)  spin  electrons  than  majority  (up)  spin.  This  is  quite  pronounced  in 
the  case  of  the  f.c.c.  Fe13  cluster,  where  there  is  a  minority  spin  excess  on 
the  central  atom  of  about  2  electrons.  In  a  sense,  the  Fe13  cluster  can  be 
reyarded  as  similar  to  an  antiferromagnet.  But  in  contrast,  the  nickel 
clusters  studied  do  not  exhibit  this  effect:  both  the  value  of  the  average 
moment  and  its  spatial  distribution  are  closer  to  that  of  the  bulk  metal  than 
for  iron. 

These  results  are  in  accord  with  the  ideas  of  spin  density  functional 
theory  in  regard  to  the  tendency  of  the  exchange  correlation  potential  to 
induce  a  spatial  separation  of  majority  and  minority  spin  densities.  This 
tendency  is  much  more  pronounced  in  iron  clusters  than  in  nickel  because  the 
atomic  moment  is  larger. 

The  ionization  potentials  of  the  Feg  and  Fe^3  clusters  have  been 
determined  by  separate,  self  consistent,  transition  state  calculations. 

In  an  independent  calculation  of  basic  importance  to  density  functional 
theory,  we  have  calculated  the  temperature  dependence  of  the  exchange 


correlation  potential  for  a  spin  polarized  electron  gas  over  a  wide  range  of 
temperatures.  Relatively  low  temperatures  are  relevant  to  solid  state 
physics,  and  in  this  situation  we  find  that  the  leading  (T^)  temperature 
dependence  of  the  exchange  and  correlation  potentials  cancel:  over  the  range 
of  temperatures  of  interest  to  solid  state  physics,  the  temperature  dependence 
of  the  exchange  correlation  potential  is  negligible  except  possibly  in 
circumstances  in  which  the  electron  density  is  exceptionally  low. 
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